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Abstract 

We argue that the analysis of Reimers et al. [ Phys. Rev. B 45, 7295 
(1992)] of their Monte Carlo data on the Heisenberg pyrochlore antiferromag- 
net, which suggests a new universality class, is not conclusive. By re-analysis 
of their data, we demonstrate asymptotic volume dependence in some ther- 
modynamic quantities, which suggests the possibility that the transition may 
be first order. 
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Reimers, Greedan and BjorgvinssonEl have recently presented results of Monte Carlo 
simulations of the pyrochlore Heisenberg antiferromagnet. An important issue, as they have 
emphasised, is to determine the order of the phase transition. This highly frustrated system 
is interesting from this point of view since Azaria et a/.i have conjectured that frustration of 
this type can easily induce a first-order transition. Unfortunately, the determination of the 
order of a transition by numerical simulations on finite systems is not always straightfoward, 
as emphasized in Ref. [ p]]. Peczak and Landau,i note that if the correlation length for a 
system with a very weak first-order transition is much larger than the system size, then 
pseudocritical behavior may be observed. Scaling of the minima near T/v in Binder's fourth- 
order energy cumulanti,V =1 — ^ < E 4 > / < E 2 > 2 , with system size L (for, e.g. an 
LxLxL lattice) is believed to be useful for the characterization of the order of a transition; 
however, a decisive conclusion based on this quantity alone can be precarious. This is due 
to the fact that one can never, by simulations of finite systems, determine if V^ in is exactly 
| or less than this value, as expected for continuous and first-order transitions, respectively. 
A value slightly less than | may indicate either a continuous or very weakly first-order 
transition. Nevertheless, it is of interest to note that the 5-state Potts model in 2D, which 
is considered to exhibit one of the weakest first-order transitions known, has a correlation 
length estimated to be on the order of 2000 lattice spacings and meaningful results have been 
obtained from simulations of lattices that were a factor of 10 smaller!!. The extrapolated 
value, V£ in ~ 0.6655, is consistent with expectations. 

In addition to V^ in , Reimers et aM considered the scaling behavior of the specific heat 
and susceptibility maxima near T N , C max (L) and Xmax(L). The pyrochlore lattice contains 
16 spins per unit cell and simulations were performed on lattice sizes L=3-10 (with L=9 
omitted). In terms of more familiar Bravais lattices with one spin per unit cell, the number 
of sites involved correspond approximately to L=7.6 — 25.2. It is certainly not clear if 
these lattices are large enough to reveal the true critical behavior of a highly frustrated 
system. Using data for all lattice sizes (except in the case of V min (L) where L=3 results 
were excluded), they found exponent ratios estimated to be a/z/=1.8(2) and 1.7(2), from 



the scaling of V m in and C max , respectively, as well as 7/V=3.5(3) from results for Xmax- In 
the case of a first order transition, these thermodynamic quantities should scale as L 3 . They 
also made the extrapolation V^ in =0. 66658(10) from their data. With these values, as well 
as other scaling results for (3 and u, they conclude that the transition is continuous and not 
within any known universality class. 

Figs. 1, 2 and 3 present an alternative scaling analysis, namely as a function of L 3 , of 
their data for \ — V m i n -, C max and Xmax, respectively. The results suggest that the asymptotic 
behavior of these quantities are in fact consistent with the transition being first order. It 
appears possible that it is only for L larger than about 6 or 7 do the simulation results 
become meaningful, as indicated by the best-fit straight lines to data points for these lattice 
sizes. The obvious exceptions are the specific heat and energy cumulant result for L=8, for 
which we have no explanation. Our best-fit estimate for V£ in from Fig. 1 is 0.66643(10), 
noteably further than from the value | than the result of Reimers et a/.S. 

From this analysis, it is clear that the hypothesis of a first order transition associated with 
the pyrochlore antiferromagnet can be taken seriously. The results of further simulations on 
larger lattices would be necessary in order to make a more definitive statement regarding 
the critical behavior of this system. It is unlikely that experimental work can resolve this 
issue. Even if the transition is in fact weakly first order, there may well be a significant 
region of temperature near T N where pseudocritical exponents can be measured, with results 
consistent with those estimated by Reimers et alK This was indeed found to be the case 
for the exponent (3. 

This work was supported by NSERC of Canada and FCAR du Quebec. We thank J. 
Reimers for useful discussions. 



3 



REFERENCES 

1 J.N. Reimers, J.E. Greedan, and M. Bjorgvinsson, Phys. Rev. B 45, 7295 (1992). 

2 P. Azaria, B. Delamotte, and T. Jolicoeur, Phys. Rev. Lett. 64, 3175 (1990); J. Appl. Phys. 
69, 6170 (1991). 

3 P. Peczak and D.P Landau, Phys. Rev. B. 39, 11932 (1989). 

4 M.S.S. Challa, D.P. Landau, and K. Binder, Phys. Rev. B 34, 1841 (1986). 
5 M. Fukugita, H. Mino, M. Okawa, and A. Ukawa, J. Phys. A 23, L561 (1990). 



4 



FIGURES 

FIG. 1. Scaling of the energy-cumulant minima with volume. The straight line is a fit to the 
data L=6-10 

FIG. 2. Scaling of the specific-heat maxima with volume. The straight line is a fit to the data 
at L=6, 7, and 10. 

FIG. 3. Scaling of the susceptibility maxima with volume. The straight line is a fit to the data 
L=7-10 
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